Spin-Dependent Ringing and Beats in a Quantum Dot System 
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We report spin-dependent quantum coherent oscillations (ringing) and beats of the total and the 
spin currents flowing through a quantum dot with Zeeman split levels. The spin dependent transport 
is calculated via nonequilibrium Green function in the transient after a bias voltage is turned on 
at i = 0. The dot is coupled to two electrodes that can be ferromagnetic or nonmagnetic. In the 
ferromagnetic case both parallel and antiparallel alignments are considered. The coherent oscillation 
and beat frequencies are controlled via the Zeeman energy Ez- In particular, for Ez = no beats 
are observed and the spin current is zero for nonmagnetic leads. In the ferromagnetic case a finite 
spin current is found for Ez = 0. The effects of temperature are also analyzed. We observe that 
with increasing temperature the ringing response and the beats tend to disappear. Additionally, 
the spin current goes to zero for nonmagnetic leads, remaining finite in the ferromagnetic case. 
The tunnel magnetoresistance (TMR) also reveals quantum coherent oscillations and beats, and it 
attains negative values for small enough temperatures and short times. 

PACS numbers: PACS number 



I. INTRODUCTION 



Spin coherent dynamics and transport in quan- 
tum dots has attracted a lot of attention due to 
its relevance to the potential new generation of spin- 
tronic devices^ (e.g., Datta-Das transistor— and storage 
devices^i^i^) , and for quantum computation and infor- 
mation processes where quantum coherence of the 
electron spin is desirable. Recent experiments demon- 
strate the possibility to coherently manipulate quantum 
states of single and double electron spins in quantum dot 
systems "'^^'"'^"'^i-'^^'-'^'^i-'^^i"'^^'"'^^ Such control can be achieved, 
for instance, via fast bias/gate voltage pulses, electron 
spin resonance (ESR) fields and coherent optical fields. 
Those experiments indicate the feasibility of using a sin- 
gle electron in quantum dot system as a quantum bit and 
reveal encouraging spin coherence lifetimes for quantum 
processing. 

In addition, a variety of interesting spin-coherent ef- 
fects have been reported during recent years. For 
instance, coherent Rabi oscillations generated by 
an ESR field were studied both theoreticallyii and 
experimentallyi^ in quantum dot systems in the Coulomb 
blockade regime. Coherent spin-dynamics was probed 
via time-resolved Faraday rotation (TRFR) technique, 
revealing coherent dynamics of net carrier spins in semi- 
conductor bulk, quantum well and quantum dotSj ^^'^^i^° 
In addition to quantum coherent oscillations, quantum 
beats of the spin magnetization due to Zeeman split- 
ting in a semiconductor quantum dot system were ob- 
served via time-resolved Faraday rotation techniquCfi^ 
and quantum beats in the TRFR signal were measured 
for magnetic fields higher than 5T in a quantum dot 
systemJ^ Quantum coherent beating in the magnetore- 
sistance of a double-barrier structure with both diluted 
magnetic semiconductor (DMS) well and contacts was 
also predictedi^i Additionally to this, coherent oscilla- 





FIG. 1: System studied: a quantum dot attached to two leads 
via tunnel barriers. Both ferromagnetic and nonmagnetic 
leads are considered. The dot level is Zeeman-split due to 
a external magnetic field. Before the bias voltage is turned 
on (t < 0) the spin split level is energetic accessible for the 
electrons in the leads, while (> Ef) is forbidden. After a 
bias voltage Vo is applied {t > 0) both levels become inside 
the conduction window and a spin-dependent tunnel current 
arises. 



tions of polarized current in magnctotransport through 
Zeeman split levels of quantum dots were investigated via 
quantum rate equationsi^ 

Quantum coherent oscillations (ringing) of the current 
fiowing through a single level quantum dot attached to 
nonmagnetic leads was predicted in a transient time scale 
after a bias voltage is turned onj ^'^i^^i^^i^^ Here we de- 
scribe transport in the same transient regime but includ- 
ing spin-dependent effects. An understanding of the spin- 
full case can better guide possible experimental measure- 
ments of the ringing response which now comes up with 
new spin-based signatures. Here we apply the nonequi- 
librium Green function (NEGF) formulation developed 
in 01 to account for (i) Zeeman splitting of the dot level 
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and/or (ii) ferromagnetic leads. These two additional 
features introduce new spin-dependent effects. For ex- 
ample. the coherent oscillations reported previously in 
j23ll24|[25 t become spin-dependent with slightly different 
frequencies for each spin component of the current. This 
spin-splitting of the frequency generates quantum coher- 
ent beats of the total current 1 = 1-^+11 and of the spin 
current I" — I-^ — I^. We also find a dynamical tunnel 
magnetoresistance (TMR) that shows coherent oscilla- 
tions, coherent beats and attains negative values in the 
transient regime for small enough temperatures. 

The paper is organized as follows: in Sec. II we de- 
scribe the system and the formulation adopted, in Sec. 
Ill we present and discuss the results and in Sec. IV we 
conclude. 

II. MODEL AND FORMULATION 

Figure [1] illustrates the system considered. It is com- 
posed of one quantum dot coupled via tunnel barriers to 
a left (L) and to a right (R) lead. Both ferromagnetic 
and nonmagnetic leads are considered. For the ferro- 
magnetic case both parallel (P) and antiparallel (AP) 
alignments are analyzed. We assume that for zero bias 
voltage {t < 0) the spin-f and spin- J, levels are above and 
below the Fermi energy Ep of the reservoirs, respectively. 
This prepared configuration gives rise to a relatively high 
transient spin current in the case of nonmagnetic leads 
as we will describe in Sec. Ill(a). When the bias voltage 
is turned on {t > 0) both levels attain resonance with 
the emitter states, thus resulting in the subsequent spin 
dynamics. 

The total Hamiltonian of the system is given hy H = 
Hl + Hr + Hd + Ht, where 



H 



L/R 



>^kaL/R{t)cl^j^/lfC]^aL/B., 



(1) 



with ekaL/R{t) = EkaL/i? + ^L/R{t) being the free- 
electron energy with a time-dependent contribution 
^L/Rit) due to the bias voltage. The operator C]^aL/R. 
(''ko-L/ij) a-iinihilates (creates) one electron with wave vec- 
tor k and spin a in the left (L) or right (R) lead. For the 
dot Hamiltonian we have 



Hr 



(2) 



where ed[t) — + ^d{t) is the dot level with a time- 
dependent term Ad{t), Ez is the Zeeman energy due to 
some external magnetic field^l and cr = + or — for spin 
t or 1, respectively. We neglect the Zeeman splitting of 
the leads E^^^'^'^ in our modeli^ The operators and 
annihilate and create, respectively, one electron with 
spin (T in the doti^ For the leads-dot coupling we have 



Ht = Y^^Vcl^d^ + V*dUkan), 



(3) 



War] 



where V is a. constant coupling parameter. We do not 
account for spin relaxation and spin dccohcrcncc in our 



model, which is reasonable for short enough times. Pre- 
vious studies have found for the relaxation and dccoher- 
ence times, Ti and T2 respectively, typically Ti 1 — 20 



ms 
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and T2 > ^s 
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The time dependent current is given by a sum of the 
currents flowing into and out of the dot, rj{t) ~ rj'™{t) + 
In the nonintcracting case and wideband limit 
they read^ 



eVl 

h 



de 



Me)lm[A,^{e,t)], (4) 



and 
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en 



i=L.,R 



(5) 



where F^ is the tunneling rate, F^ — 2n\V\'^p^l, with 
being the constant density of states for spin a in lead 
T]{= L, R). The ferromagnetism of the leads is accounted 
for via the model F^ = Fo(l -h ap) and F^ = Fo(l ± ap), 
where Fg is the leads-dot coupling strength and p gives 
the polarization degree of the leads. The signs +/— 
in F^ apply for parallel and antiparallel configurations, 
respectively!^ The function /,,(e) is the Fermi distribu- 
tion function for lead rj. For a bias voltage of the kind 
Al = 0, Ar = -Vo0{t) and Ad = Aii/2 the function 
A<j,,(e,t) is given by2^ 



A^^(e,t > 0) 
1 - e 



„i[e-eo-o--Ez/2-Ad+A^+i(r^+r^)/2]t/R 



e-eo-(7Sz/2 + ^(F^+F^)/2 

i[£-£o-crBz/2-Ad + A„+j(r^+r«)/2]t/'i 



e-eo- <7Ez/2 - A^ + A^ + i{TL + F^)/2 ' 



(6) 



Using this expression inside Eqs. ^ and ^ we find the 
results presented in the next section. 

For comparison we apply the master equation (ME) 
technique to calculate the current and spin-current. 
While the NEGF is vaHd for both fesT > Fq and ksT <C 
Fq, the ME approach is accurate only for ksT ^ Fg 
(sequential-tunneling limit). So we will restrict the com- 
parison only for relatively high temperatures, which cor- 
respond to the range of mutual validity of the two ap- 
proaches. The current expression in the ME formulation 
is given by22 

where Pqi Pa and P2 are the probabilities to have no 
electron, one electron with spin cr and two electrons, re- 
spectively, in the dot. The Fermi functions are /,,cr = 
{l+exp[(e,-A^)/(fcsT)]}-i and = {l+cxp[(e,+;7- 
A,,)/(fcBT)]}^^ where Ecr = ed+crEz /2. Defining the vec- 
tor of the occupation probabilities P = (Pq, P|, P|, P2)^ 
we can write the master equation as 



P = M P, 

where the transition matrix is given by M 
with 



(8) 



R 
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M" = 



v f t 



-r;'(i 



1 - M 

-M-^Jfvi 







-r;%T-rj(i-/,,] 







rT'(i-/„T) 
r;%T + n'^-(q + rp/ 



(9) 



If we take fucr = fna = we obtain a matrix form for M 
similar to the one presented in [3^ . Even though Eqs. 
([7])- ([9]) can account for Coulomb interaction in the se- 
quential limit, we are interested in comparing the results 
obtained from both the ME and the NEGF (which is non- 
interacting in the present formulation), thus we simply 
assume [/ = in the above ME expressions. Some effects 
of Coulomb interaction in a transient response (sequen- 
tial limit) can be found, for instance, in [39j. 

III. RESULTS 
A. Nonmagnetic Leads 

Figure [5] shows the time evolution of the total cur- 
rent / = /| + /| and the spin current I'* = /| — /| in 
the emitter lead;^ obtained via NEGF. Different Zee- 
man energies are considered. For Ez = there is no spin 
current [Fig. [IJb)] and the current [Fig. [DJa)] presents 
the typical coherent oscillations that arise after a bias 
voltage is turned 0^, 23,24,25,26 ^pj^g oscillatory frequency 
is given by fnoc = \Ef — Cd\ (or equivalently the period 
Tc = 2iTh/\EF — Cdj) and the damping is due to the 
leads-dot coupling. We may note that for the parame- 
ters adopted in Fig. [2ja) we have Tc « 0.06ft/ro4i In 
the presence of Zeeman splitting the frequency ojc be- 
comes spin dependent, with hujc ~ \Ef — Cd — aEz/2\, 
so the spin components of the current oscillate with dis- 
tinct frequencies. In the case of relatively close frequen- 
cies, quantum beats of / and are seen, with beating 
frequency given by h\Lul — u;^| = Ez [Figs. [l{c)-(f)]. 

The spin current observed in the presence of Zeeman 
energy [Figs. [2l^d) and (f)] comes from the initial config- 
uration {t < 0) of the system. For t < the occupations 
are w 1 and n-f w 0. When the bias voltage is turned 
on at i = the spin j electrons in the emitter lead can 
start to flow through the dot, while the spin J, ones have 
to wait till the initial J, electron in the dot leaves to the 
collector lead (Pauli blockade). Since this tunnel event 
takes a typical time of H/Tq a spin current is expected 
in this time range. In the case of initially both 7i| « 
and 7i| « 0, which can be achieved by taking eo high 
enough in order to forbid thermal occupation of the spin 
split levels, an oscillatory spin current around zero is ob- 
served. 

In Fig. ([3]) we show /| and /| separately. For Ez = 
the currents /| and /| coincide, while for Ez = 10 
and Ez = 20 they differ, thus generating spin current. 




0.0-1 . . . . . o.oJ : ^ ^ : 

0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

time time 



FIG. 2: Total current 7| -|- 7| (left panels) and spin current 
/| — 1 1 (right panels) against time for different Zeeman ener- 
gies Ez- For Ez = the current shows coherent oscillations 
(ringing) with a period Tc — 2nh/\EF — edj. No spin current 
is observed for Ez = 0. In contrast, for Ez = 10 and Ez = 20 
a spin current arises and quantum coherent beats are seen in 
both current and spin current. Units: eFo/h for the currents, 
h/To for the time and To for the energies. Parameters: eo — 0, 
T/o = 200, ksT = 0.1, p^O. 



Eventually for high enough times {t ^ 1) both /| and 
become close to each other, thus resulting in 0. The 
difference in frequency for each spin component is clearly 
seen in Figs. [3ljb)-(c). In particular as Ez increases, we 
observe a blue shift for the frequency oj^ and a red shift 
for ojI. 

In Fig. (|3]) we analyze the effects of temperature on 
the coherent oscillations. We find that as ksT increases 
the oscillatory behavior tends to disappear in both the 
current and the spin-current. This is expected since with 
increasing T the coherence is washed out. For comparison 
we show in Fig. me)-(f) the total and the spin currents 
calculated via the master equation technique (dots). We 
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FIG. 3: Spin resolved currents (denoted by up and down ar- 
rows close to the curves) against time for three different Zee- 
man energies. For Ez = 0, 7^ and ly coincide. As Ez in- 
creases we find a red-shift for the frequency of and a blue- 
shift for the frequency of Jj. This frequency splitting results 
in the quantum beats seen in Fig. [5] Units and parameters 
as in Fig. [2l 



find close agreement between the nonequilibrium Green 
function and the master equation results, except by the 
strongly suppressed coherent oscillations still seen in the 
NEGF calculation for ksT — 10. As the temperature 
increases these coherent oscillations (ringing) disappear 
and both NEGF and ME approaches give equal results. 
Additionally to this, we note that the spin current is sup- 
pressed with increasing ksT. For even higher tempera- 
tures we find zero spin current for all times, which is re- 
lated to the fact that n^{t < 0) approaches ni{t < 0) with 
increasing ksT-thns resulting in a less effective tran- 
sient Pauli blockade. When the leads arc ferromagnetic, 
though, the source of spin current is not only the initial 
configuration of the system but also the leads itself, so 
J'* remains nonzero even for high temperatures, as we 
describe next. 



B. Ferromagnetic Leads 

Figure ([5]) shows the current and the spin current 
against time in the presence of ferromagnetic leads for 
differing temperatures and Ez = 10. Both P and AP 
configurations are analyzed. Here we have enlarged the 
time scale compared to Fig. |4] because some contrasting 
features between P and AP alignments are better seen 
after t = 1 (e.g., Ip > Iap)- For t < 1 the standard 
inequality Ip > Iap is only true when we approach the 
sequential limit for fc^T = 10. Curiously, for fc^T = 0.1 
and ksT = 1 wc find Iap slightly above Ip for t < 1. 
This unusual behavior can be qualitatively understood 
in terms of tunnel rates and the initial (prepared) spin- 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

time time 

FIG. 4: Total current (left panels) and spin current (right 
panels) against time for differing temperatures. Quantum 
coherent oscillations and beats are seen for fcsT = 0.1 and 
ksT = 1. As the temperature increases the oscillations vanish 
(see ksT — 10). Additionally, the spin current is significantly 
suppressed for ksT = 10. For comparison in panels (e) and 
(f) we show the results obtained via master equation (dots). 
Units as in Fig. [g] Parameters: eo = 0, Vo = 200, Ez = 10, 
p = 0. 



population of the dot. Initially {t < 0) the dot is occupied 
by one spin J, electron. When the bias voltage is applied 
(t > 0) this electron can leave the dot to the collector 
lead with a tunnel rate given by T^. Since F^^ is big- 
ger for AP than for P configuration, the initial spin | 
electron leaves the dot faster in the AP alignment. This 
turns the Pauli blockade less effective in this configura- 
tion, thus allowing the emitter current Iap to be slightly 
greater than Ip at initial times. 

When ksT exceeds the coupling strength Tq, the co- 
herent oscillations of both the current and the spin cur- 
rent are strongly suppressed and the standard inequality 
Ip > Iap is recovered even for short times after t = 0. 
Here we also see an agreement between the results ob- 
tained via NEGF and ME, except by the coherent os- 
cillations residually present in the NEGF calculation for 
ksT = 10. As the temperature increases even further 
these residual oscillations of the current and the spin- 
current are washed out and both NEGF and ME ap- 
proaches converge to the same results. Finally, we ob- 
serve that the spin current does not tend to zero with in- 
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FIG. 5: (Color online) Total current / = /| + /| (left pan- 
els) and spin current I" — — J| (right panels) against 
time for differing temperatures in the case of ferromagnetic 
leads. Both parallel (P) and antiparallel (AP) configurations 
are shown. Coherent oscillations and beats are seen in both 
alignments. The standard inequality Ip > Iap (magnetore- 
sistance effect) is seen for t > 1 and kpT = 0.1 or ksT = 1. 
For times t < 1 unusual behaviors like Ip < Iap can be found 
for these temperatures. In contrast, for kpT — 10 almost no 
oscillations are seen and Ip > Iap is recovered for all times. 
For comparison in panels (e)-(f) we show results obtained via 
master equation technique (dots). We also note that the spin 
current goes to zero as the time evolves in the AP alignment 
while remains finite in the P case even for kpT = 10. Units 
as in Fig. [1 Parameters: eo = 0, Vb = 200, Ez = 10, p = 0.4. 



creasing temperature as we see in the nonmagnetic case 
[compare Figs. IUf) andj^f)]. The suppression of I'* is 
observed only with increasing time in the AP ahgnmcnt, 
remaining finite for the stationary limit in the P config- 
uration. 

In Fig. ^ we show the time-dependent TMR, de- 
fined as TMR=(/p — Iap) /Iap- Coherent oscillations 
and quantum beats are seen for both fc^T = 0.1 and 1 
and is only residually observed for fc^T = 10. Interest- 
ingly, the TMR attains negative values for short times 
and small temperatures4^ In panel (c) we compare the 
TMR obtained via NEGF (solid line) and ME (dots) with 
essentially no contrast between them. 

Finally in Fig. [7] we show the spin resolved cur- 
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FIG. 6: Tunnel magnetoresistance (TMR) against time for in- 
creasing kpT. For kpT = 0.1 and 1 the TMR shows a ringing 
response, quantum beats and attains negative values before 
it increases to a stationary value close to 20%. In contrast, 
for kpT = 10 the TMR is always positive with almost no os- 
cillations. The dots show the TMR obtained via the master 
equation technique. Units as in Fig. [21 Parameters: eo — 0, 
Vo = 200, Ez = 10, p = 0.4. 



rents and the spin current in the ferromagnetic case for 
Ez ~ 0. Here 7^ [in contrast to the nonmagnetic 
case, Fig. [H^b)] and no beats are seen since /| and 
have the same frequency. Additionally, the oscillation 
amplitudes and the decaying rates can differ for each 
spin component. In the parallel configuration we observe 
that /| starts oscillating with a higher amplitude than 
/| but its oscillations are faster suppressed. This is due 
to the inequahty -f r|*^ > Fj' -I- r|*. The interplay 
between oscillation amplitudes and decaying rates gives 
rise to a node in the spin current seen around t = 0.5 
[Fig. Eljb)]. For long enough times both /| and /| at- 
tain their respective stationary values with /| > /| [Fig. 
[7{a)]. In the AP alignment /| starts oscillating with a 
higher amplitude than /| but, in contrast to the P case, 
the decaying rate is the same for both spin component 
(F^' 4-Fj^ = F| +rf )• Consequently no node is observed 
in the spin current, which simply oscillates decaying to 
zero. 



IV. CONCLUSION 

We have calculated spin-dependent transport through 
a Zeeman split quantum dot level after a bias voltage 
is turned on. The dot is coupled to two leads (source 
and drain) that can be ferromagnetic or nonmagnetic. 
Quantum coherent oscillations (ringing) and beats of the 
current and of the spin current are found in a transient 
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time scale. The frequency of the beats can be tuned via 
the Zeeman sphtting energy Ez- We also report a spin 
current that arises when Ez ^ and fc^T < Fq for 
nonmagnetic leads. For the ferromagnetic case the spin 
current can be seen even for Ez = and high tempera- 
tures (/cfiT > Fq). In particular /" goes to zero as the 
time evolves in the AP alignment and remains finite in 
the stationary limit for the P configuration. The tunnel 
magnetoresistance was also analyzed. We found quan- 
tum coherent beats and negative values of the TMR due 
to the transient dynamics of incoming and outgoing spin 
polarized electrons in the quantum dot. The negative 
vales of the TMR vanish as the temperature increases, 
which is confirmed by master equation calculations. 



FIG. 7: (Color online) Spin resolved currents (denoted by up 
and down arrows closed to the curves) and spin current for 
both P and AP alignments and Ez ~ 0. Both /| and 
oscillate with the same frequency due to Ez = 0, but in con- 
trast to the nonmagnetic case [Fig. EJa)] here the amplitude 
and decaying rate (F^ -I- F^) can differ for each spin compo- 
nent. In particular in the P alignment the interplay between 
the amplitudes and the decaying rates gives rise to a node in 
the spin current. Units as in Fig. [5] Parameters: eo = 0, 
Vo = 200, Ez = 0, ksT = 0.1, p = 0.4. 
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